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iHj! Abstract 

We discuss various aspects of the inner structure formation in viri- 

■ alized dark matter (DM) halos that form as primordial density in- 
00 ! homogeneities evolve in the cosmological standard model. The main 

■ focus is on the study of central cusps/cores and of the profiles of DM 
(<~^ ■ halo rotation curves, problems that reveal disagreements among the 
Q,^ . theory, numerical simulations, and observations. A method that was 
Q I developed by the authors to describe equilibrium DM systems is pre- 

■ sented, which allows investigating these complex nonlinear structures 
. . , analytically and relating density distribution profiles within a halo 

, ^ I both to the parameters of the initial small-scale inhomogeneity field 

^ ■ and to the nonlinear relaxation characteristics of gravitationally com- 

pressed matter. It is shown that cosmological random motions of mat- 
ter 'heat up' the DM particles in collapsing halos, suppressing cusp- 
like density profiles within developing halos, facilitating the formation 
of DM cores in galaxies, and providing an explanation for the differ- 
ence between observed and simulated galactic rotation curves. The 
analytic conclusions obtained within this approach can be confirmed 
by the N-body model simulation once improved spatial resolution is 
achieved for central halo regions. 
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1 Indroduction 

This review is the second of two (see [I]) devoted to the formation of the 
large-scale structure of the Universe and problems of virialized dark mat- 
ter (DM) halos. We do not attempt to discuss all related problems of the 
physical cosmology and focus in fact on DM evolution in the cosmological 
standard model (CSM). The discussion is based on the original papers of 
the authors and follows the chapters of book [2j. Here, we present a more 
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detailed discussion of these problems using new observational data and a 
comparison with the theory without invoking detailed analytic calculations 
(theoretical results are described in the appendices; see |2] for the proofs). 

2 The essence and history of the problem 

Explaining the internal structure formation of virialized DM halos, which are 
studied both observationally and using analytic and numerical simulations, 
is one of the key cosmological issues. 

Despite cosmologists' significant efforts, many aspects of halo formation 
remain unclear. In particular, it is unclear whether the universal density pro- 
file [3]-[S] (hereafter, the Navarro- Frenk- White (NFW) profile), discovered in 
a wide range of masses and sizes of DM halo formation models, is truly uni- 
versal. The NFW profile has been found to correspond well to the observed 
density profiles in galaxy clusters [6], but at the same time it deviates from 
the observed galactic profiles in many respects (see [ZIIH])- This means that 
many physical factors controlling galaxy formation are not fully taken into 
account in numerical models. 

Numerical calculations of the N-body problem, which have been widely 
used for many years, are presently a powerful tool to study the formation of 
the nonlinear structure of the Universe. But numerical models have many 
restrictions, which stimulates the development of analytic and approximate 
methods that can reveal and qualitatively characterize the effect of the main 
physical processes responsible for halo formation. One of the key processes 
that are difficult to model numerically is the relaxation of matter contracted 
in a protohalo. To analyze this process, the cosmological small-scale matter 
velocity and density fields should be taken into account. 

In the framework of the widely used spherically symmetric model of the 
collapse, direct relations between the properties of the initial density per- 
turbations and virialized halos are given by the Press-Schechter relations or 
their extensions (see [9]- [2^). In this approach, the collapse is described as 
a consecutive fall of spherical shells, with their subsequent relaxation and re- 
distribution inside the halo. Each spherical layer is characterized by the mass 
and turn-around radius (the moment of 'decoupling' from the Hubble expan- 
sion), which can be related to the spectrum of initial density perturbations, 
for example, using the smoothing of the density perturbation field by filters 
of different scales (see \i7\ ) . This approach can be generalized to include the 
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angular momentum of the shells, external tidal forces, an external density 
distribution, etc. [2T]-|2B]. Another description of halo formation,which is not 
based on the quasispherical collapse but uses the Zeldovich approximation, 
is proposed in |27] . 

It is well known that the formation of any dense object means the re- 
laxation of contracted matter. In the spherical collapse model, this is the 
so-called violent relaxation process, which provides the energy and mass re- 
distribution of the compressed matter, the removal of the excess energy from 
the system, and the halo contraction. The simplest example of the violent 
relaxation is given by the collapse of a homogeneous nonrotating DM cloud 
of a finite size (see [2H]-[3n])- The angular momentum does not stop the 
violent relaxation [2T]-[23], but does weaken it (see, for example, [SBf IST]). 

In contrast to such methodically transparent results, numerical calcula- 
tions show that to describe the halo formation process, it is insufficient to con- 
sider the fall of only diffusive matter (spherically symmetric or anisotropic). 
In numerical models, the formation of a halo is accompanied by the formation 
of random systems of low-mass high-density subhalos and their subsequent 
coalescence in the main protohalo. This process of hierarchical merger pro- 
vides the mass transfer inside the main halo, in which small subhalos are 
'stripped' by the tidal interaction and dynamical friction with the surround- 
ing matter (see [32]- [35]). Both processes - the fall of diffusive matter and 
the merger of satellite subhalos - play similar roles in the halo formation. 

The results of numerical calculations suggest that the merger of early 
formed high-density subhalos can stimulate the formation and later increase 
in the central density in the main halo, and the violent relaxation here has 
a smaller effect. For example, even one merger of subhalos with comparable 
masses can strongly change the internal structure of the whole halo [5B] . 
The merged subhalos supply the common halo with their kinetic energies and 
angular momenta, and hence the destruction and stripping of subhalos inside 
the main halo become the principal constituents of the relaxation process. 
The interaction of subhalos with the main halo, other subhalos, and diffusive 
matter redistributes the energy in the main halo and changes its internal 
structure. The efficiency of relaxation in numerical calculations is comparable 
with estimates obtained in the simplest spherically symmetric collapse models 
(see [an Ian [38]). 

The relaxation is needed to provide the energy removal from the main 
halo, but at the same time it leads to the chaotization of the collected ma- 
terial and largely determines the density profile unification in the relaxed 
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halos. Nevertheless, in both the hierarchical and spherical halo formation 
models, the merger of subhalos (or spherical shells) is not purely chaotic and 
is controlled by the initial density and velocity distributions. During the re- 
laxation, the initial conditions do not smear out completely and are partially 
conserved. As we show below, both halo mass distribution and their internal 
structure depend on the initial perturbation spectrum. 

It is well known that initial conditions can dramatically affect the evo- 
lution of dynamical systems. The difference between the cold and hot dark 
matter cosmological scenarios provides an example. In the context of halo 
formation, the possible influence of thermal velocities and small-scale ini- 
tial perturbations were discussed in [221 EZl ESI SO]- Here, we significantly 
improve the description of this effect by using statistical characteristics of 
spatial density and velocity distributions inside a collapsed DM cloud |¥T] . 
This approach allows estimating the background entropy of the relaxed mat- 
ter, which is related to the initial small-scale velocity perturbations, and its 
dependence on the mass of halos formed. The effect of factors such as the an- 
gular momentum of collapsing matter and nonlinear destruction of low-mass 
subhalos is not considered here. Due to the random character, these factors 
can only increase the background entropy and its deviation from the mean 
value. Our calculations yield the minimum level of the background entropy, 
which, nevertheless, can change the internal structure of the emerging halos 
and significantly suppress the central DM cusp formation inside them. 

The use of the entropy approach in describing relaxed halos allows the 
determination of the total entropy (initial entropy and that acquired during 
nonlinear relaxation) of DM particles in the halo. In contrast to the kinetic 
approach, our method does not appeal to the halo particle distribution and 
does not describe the relaxation history. But it correctly reproduces any dis- 
tribution of matter inside a relaxed spherical halo and distinguishes between 
the properties of this distribution that are due to the initial conditions and 
those that are related to the relaxation process itself. Such a calculation is 
possible due to the property of entropy as the most stable characteristic of 
matter, changing only during irreversible processes. 

For both spherical and hierarchical halo formation models, it is possible to 
estimate the entropy growth during matter relaxation [42j . For the spherical 
collapse, this estimate can be done analytically (see [28]). The hierarchical 
relaxation process is more complicated and is not described analytically, and 
we therefore use the results of numerical calculations to estimate the entropy 
growth there. 
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In the observed gravitationally bound objects, the spatial distributions 
of DM, baryons, and luminous matter are different. In the DM-dominated 
Universe, it is reasonable to regard the DM halo formation problem as the 
first approximation to a more complicated problem of the formation of grav- 
itationally bound systems. Here, we consider only the DM dynamics in the 
CSM framework, ignoring the effect of dark baryons and luminous matter. 
Nevertheless, in some cases where the effect of the baryon component on the 
internal structure of a galaxy can be significant (see, e.g., [43]), this issue 
requires further studies. 

3 Dark matter halo 

We assume DM to consist of nonrelativistic massive particles that interact 
with each other and other particles only gravitationally. This matter is called 
'dark' because it is invisible - it has a nonbaryonic nature and does not 
interact with light Q However, DM can be studied using dynamical methods, 
because it is clumped and produces spatial gradients of the gravitational 
potential that influence the motion of visible bodies (galaxies, stars, gas), the 
state of baryons (hot gas), and light ray deflection (gravitational lensing). 

The mean DM density in the Universe is five times higher than the density 
of cosmological baryons, which is why nonrelativistic dark particles control 
the process of gravitational clustering and of the evolution of the inhomo- 
geneous part of the gravitational potential of the Universe. The density 
contrast of a weakly inhomogeneous initial spatial DM distribution increases 
with time. Because this matter is cold, the pressure gradients there are small 
and cannot prevent the development of gravitational instability. 

In regions with enhanced density, the Friedman expansion rate slows 
down, and at some instant the DM expansion stops and turns into col- 
lapse. During dynamical contraction and subsequent oscillations of matter 
flows, processes of violent collisionless gravitational relaxation occur: parti- 
cles move in a variable gravitational potential localized in space, such that 
some particles (~ 10%) are expelled from the system and carry away the 
excess positive energy. The remaining energy of the clump is redistributed 
inside the system, and a gravitationally bound DM object called a halo is 
formed. 

^Matter can be weakly visible if its particles weakly interact with light, for example, if 
they annihilate to emit photons in the regions with high DM density in galactic centers. 
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The observed halos (galaxies, their groups, and clusters) are relaxed sys- 
tems of particles that are gravitationally bound in all three spatial directions. 
The upper limits of the halo sizes and masses are a few Mpc and lO^^M©. 
We recall that the mass of matter in the Universe in a sphere with a radius 
i? = 10 Mpc is Mio ~ 1.6 • 10^^ Mo [see Section 8, formula (21)]. Here and 
below in numerical estimates, we assume the CSM parameters 

Ho ~ 70 km s"^ Mpc"^ , (1) 

where Pm is the cosmological DM density, Pc is the critical density in the 
Universe, is the cosmological density of dark energy expressed in units 
of the critical density, Hq is the Hubble constant, and G is the gravitational 
constant. 

We consider the simplest equilibrium conditions of self-gravitating DM 
systems. The distribution function of non-relativistic particles in a spheri- 
cally symmetric halo depends on the radial distance to the halo center r and 
moduli of the radial and transverse particle momenta. In spherical coordi- 
nates, the stress-energy tensor of matter takes the diagonal form 

T; = p diag (1, , (2) 

where p — p{r) is the halo density profile, and (Tr{r) and (Tt{r) are the radial 
and transverse velocity dispersions of particles. In the Newtonian limit, the 
identity T^.^ = implies the hydrostatic equilibrium equation for coUisionless 
particles: 

Id. 2 ^ 2 2 GM{r) 
pr^ ar r ar 

where the gravitational potential defined up to an additive constant, and 
the mass of the system M, which depends on r, are related to density as 

A* = f r'^^ = 47rGp , Mir) = in f pr^dr . (4) 

dr \ dr J Jq 

Inside the halo, DM particles move along different orbits, from radial to 
circular ones. Each orbit is characterized by the energy and angular momen- 
tum vector, which are integrals of motion in a static spherically symmetric 
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field. The density distribution at the halo center can be conveniently char- 
acterized by a power-law function, 

p oc r-" (5) 

with the exponent (slope) a = const G (0; 2.5). The particle distribution 
function depends on the evolutionary history of the halo. Numerical ex- 
periments and observations show that for most of the spherical halos, the 
particle velocity distribution is nearly isotropic, although deviations from 
isotropy can sometimes be as high as 20-30% (see, for example, [H]). 

Large-mass halos do not have enough time to form in the Universe. For 
example, systems collapsed in one or two dimensions, which are called walls 
and filaments, are typical nonlinear elements of the large-scale structure of 
the Universe. On average in the Universe, the density contrast on scales 
exceeding dozen Megaparsecs remains smaller than unity {\6p/p\ < 1), and 
we can consider only regions with an increased or decreased (in comparison 
with the mean cosmological) density of matter. 



4 Entropy of particles in the halo 

We consider an isotropic distribution of DM particles with ar = Ct = a 
in more detail. The halo equilibrium is determined by the balance of the 
effective pressure gradient of nonrelativistic matter, 

p = nT = pa'^ (6) 

and the gravitation produced by the total mass M = M{r) [see (3)], 

dp_ ^ GMjr) 
pdr ' 

where radial functions n = p/m and T are the density and effective temper- 
ature of particles with the mass m = const. By virtue of the equivalence 
principle, m cannot be determined from gravitational equations because par- 
ticles are moving along geodesies irrespective of their masses. The measurable 
variables are the density p and the particle velocity dispersion a. 

From the hydrostatic equilibrium equation, it is possible to find the radial 
matter density profile in the halo using the known velocity dispersion dis- 
tribution cr(r); vice versa, using the density distribution, we can reconstruct 
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the velocity dispersion law. For the adiabatic density distribution, a oc n^/^, 
where the coefficient of proportionality depends on the entropy distribution. 
In analogy with an ideal gas, we introduce the entropy function E = E{r) of 
the virialized halo as |2l E] 



e = -M'^V"«J^ = ;J^. (8) 



where rrip is the mass of a proton to which the halo particle mass m is 
normalized. 

The function E{r) is the measure of the total entropy of halo particles 
acquired during the full halo formation historyj^l It mainly includes the back- 
ground entropy determined by the initial small-scale flows and matter inho- 
mogeneities in the protohalo, and the acquired entropy generated during 
collisionless and hierarchical relaxation of matter at the nonlinear stage of 
the halo formation. 



5 Isothermal sphere 

Modeling the internal structure of observed halos frequently involves the 
isothermal sphere approximation, where the velocity dispersion can be as- 
sumed constant and independent of the radius: 

a{r) = ctq = const , kq = — 5- = — — — = const . (9) 

(To T 

This approximation is in good agreement with observational data in a re- 
stricted range of scales (Fig. 1). 

After multiplying Eqn (7) by and differentiating with respect to r, we 
obtain the equation for the density profile in an isothermal sphere: 

d / r'^dp 
pr^dr \p dr 

^We stress that E is a function of entropy and not the entropy itself. Here, we do 
not consider general issues of the applicability of the notion of entropy (coarse-grained 
entropy, see footnote 3) to collisionless particles, by defining the function E{r) only for 
stationary equilibrium systems with isotropic particle velocity distributions. We recall 
that in definition (8), cr^ = , where v is the one-dimensional peculiar velocity of dark 
matter particles at a point r (the bar means averaging over the velocity space). 
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Figure 1: Rotation curves for 2155 galaxies (from [15]). The curves, which 
are divided into several groups in accordance with the galaxy luminosities, 
are characterized by the absolute stellar magnitude in the I-filter (shown to 
the right). i?opt is the 'optical' radius of galaxies, inside which 83% of the 
integral emission flux is contained. 



This nonlinear equation contains an attractor - a particular solution that is 
an attracting separatrix for dynamical trajectories (10). The general solution 
has the form 



p(r) 



Po 



1 + rV(2r2) ' 
2 



r < To = 
r > ro (attractor) 



(11) 



where the parameter po determines the central density of the halo core and 
ro is the characteristic size of the core. Irrespective of the value of Pq, the 
density profile and mass of the halo at r > tq are determined only by the 
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velocity dispersion of particles: 



M = ^ ^ 7 . 10- ( / Y ( Me . (12) 
kq V300 kmc-i/ V200 kpc/ " ^ ^ 

For the characteristic values of a and r, we obtain the typical mass of non- 
linear dark matter halos in the observed Universe - that of galaxy groups. 
We recall that the total mass of the Local Group, which includes the Milky 
Way and the Andromeda Nebula, is 2 ■ 10^ Mq. 

We have derived Eqns (3) and (7) for DM by assuming its dominance; 
however, these equations can be used to describe hydrostatic equilibrium of 
any matter (DM, stars, gas) in the gravitational field of the total mass M(r). 
For this, the corresponding p and p should be substituted in the left-hand 
side of Eqns (3) and (7). For example, we can rewrite the halo mass in terms 
of the effective gas temperature Tgas in the equilibrium: 

M = ^ 7 ■ 10- ( ( ^--] Mo , (13) 

G/igas VlkeVy V200kpcy °' ^ ' 

where //gas — fnp = 1 GeV is the molecular weight of the gas. Such a hot gas, 
which has been present in galaxy clusters for several billion years, is observed 
by X-ray telescopes, which allows reconstructing the gravitational potential 
distribution and estimating the total mass of clusters (see, e.g., [16]). 

We return to solution (11). The core density p — po does not relate to the 
gravity of matter: it depends only on the initial cosmological conditions that 
determined the halo formation history. The size of the core is determined 
by the product of the particle velocity and the dynamical time of the central 
density: 

'"n — — ; . 

Here, tq is the radius at which the matter self-gravity becomes significant. 
At r > To, the dark matter gravity restructures the inner halo such that the 
rotation velocities of particles cease to depend on the radius: 



Vrot{> ro) = y = V2ao = const . (14) 

Such velocity distributions are realized in many gravitationally bound 
cosmological objects. For example, flat rotation curves are observed in spiral 
galaxies (Fig. 2): circular velocities of stars and gas first increase as the 
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Figure 2: The rotation curve of galaxy NGC 6503 (data from [57])- The dark 
matter halo contribution is shown by the dashed-dotted line. 

radius increases and then become constant or change very slowly. The DM 
mass increasingly contributes at large distances and maintains the flat part 
of the rotation curve over about ten dynamical scales. 

Dwarf galaxies with low surface brightness (Fig. 3) provide another ex- 
ample. There, DM dominates starting from small radii, and a linear increase 
in the circular rotation velocity is observed, suggesting a small variation of 
the DM density in the central core [see (18)]. We see that the initial increase 
of rotation curves also corresponds to isothermal sphere (11): 



The fact that the inner circular velocities are smaller than the particle veloc- 
ities in the halo (frot < o"o) suggests that the halo particles freely, by inertia, 
move over the region r ~ ro with a constant velocity and do not feel gravity. 
To within the observational accuracy, isothermal sphere distribution (11) 




(15) 
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Figure 3: The circular rotation velocity as a function of radius in low- 
brightness surface galaxies (a) LSB F571-8 and (b) LSB F583-1. The solid 
curve shows the Burkert approximation, the dashed curve shows the NFW 
approximation, and the dashed-dotted line shows the Moore profile (data 
from [is]). 



can be approximated by the simple formula 

P - T— ^ , X = —. 16 

The flat part of the circular velocity, corresponding to the linear mass increase 
M (X r [see (14)], is observed within a finite dynamic range. The density 
distribution in such a system can be conveniently described by the profile 

^ (1 + a;2) (r + r,) ^ ^ 

where = const is a characteristic cut-off parameter. In the central part 
at r < rg, the profile is close to that of the isothermal sphere in (16); but 
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at large distances, the density decreases oc and circular velocities slowly 
decrease. The circular velocity reaches a maximum in the region with the 
local density slope exponent is a ~ 2. 

6 Internal halo structure 

The radial dependence of the circular velocity is related to the density dis- 
tribution inside the halo. Figure 3 demonstrates rotation curves of dwarf 
galaxies that can be approximated well by phenomenological profile (17) 
with Tg = To, which is referred to as the Burkert approximation 09]: 



Similar examples of galactic cores with slowly changing density are typical 
for small galaxies. So far they have not been reproduced in numerical calcu- 
lations of structure formation. 

The nonlinear relaxation of a contracting cloud leads to a diverging cen- 
tral density in the forming halo. Analytic models of the collapse [28| [29] 
consider the gravitational contraction of a spherical cloud at rest with a 
smooth density distribution. After the relaxation, a singular halo density 
profile is formed with the slope exponent a ~ 1.6 — 2. This result, which 
is also confirmed by numerical calculations [50]- [52], is closely related to the 
choice of the initial state of the contracting mass, namely, to the absence of 
random small-scale velocities and an almost spherically symmetric density 
distribution. This example very well illustrates the transformation of the ki- 
netic energy of the contracting cloud into thermal energy of the halo during 
violent relaxation. 

More realistic models of the collapse with subsequent DM relaxation ap- 
pear in numerical modeling of the evolution of a system of gravitationally 
interacting N bodies (A^ ~ 10^*^). Here, random initial particle velocities, the 
anisotropy of the contraction, and the successive merger of the main halo 
and its satellites are taken into account. The mean density profile of halos 
formed in these calculations is described by the NFW approximation [3]- [5], 
which depends on a single parameter and has the form 
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(18) 



p oc 




1 



(19) 



p oc 



r (r + rg) 
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Figure 4: Dark matter mass distribution in galaxy clusters according to the 
data in [6] . 

At r < Vg, the dependence of the mean density profile on r tends asymptoti- 
cally to a cusp-like power law with the slope a = 1. 

The NFW profile obtained in numerical A^-body models well describes the 
density distribution in galaxy clusters (Fig. 4, 5). It should be remembered, 
however, that this profile is obtained by averaging density distributions of 
many halos, and deviations of individual profiles from the mean one can 
exceed 20% (Fig. 6). We also note that technical limitations restrict the 
dynamical range of model scales, and the properties of halos with low and 
moderate masses are not properly reproduced in the calculations. Numerical 
constraints also restrict the size of the central region resolved in numerical 
models. 

In numerical models, the DM halo central density tends to infinity, but 
the mass remains finite and tidal forces are insignificant. Such an internal 
structure was called the halo central cusp. Thus, singular cusps are formed 
in numerical models, but these structures are not observed in the Universe. 

Figure 7 demonstrates the rapid density increase in a DM halo that ap- 
pears in numerical models. The comparison with observed galactic profiles is 
shown in Fig. 8. The question of whether DM cusps exist in reality is impor- 
tant not only for understanding galaxy formation; it also can be connected 
with the physics of DM particles. For example, one of the methods to search 
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Figure 5: The curves shown in Fig. 4 in variables normahzed to the mass 
and radius of the sphere inside which the mean density is 200 times as high 
as the matter density in the Universe. The sohd hne shows the NFW ap- 
proximation. 

for possible DM particle annihilation assumes that central cusps do exist in 
galactic centers. 

To summarize, the comparison of the observed galaxies with numerical 
models leads to the following conclusions. 

(1) Numerical models reproduce the formation of massive objects, and 
the density profile of galaxy clusters is in good agreement with calculations. 

(2) Burkert profile (18) is mainly observed in low- and moderate-mass 
galaxies dominated by DM. The formation of such galaxies is related to the 
initial perturbations of comparatively moderate scales. In the central regions 
of more massive spiral galaxies, the DM and baryon densities are comparable 
and it is difficult to separate them. 

(3) These issues are closely related to the number count of our Galaxy 
satellites. Until recently, the number of the observed satellites was about two 
times less than predicted by numerical models |44j. However, many very faint 
satellites have been discovered recently, and their number even exceeds the 
prediction of numerical models (see [56]- [60]). This shows that the properties 
of low-mass objects need to be studied further. 
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Figure 6: The density profiles of tliree lialos witli masses 10^^ — 10^^ obtained 
in numerical models ^53j- The solid curve in the center shows the NFW 
profile. The short line in the upper part of the figure corresponds to the 
slope a = 2. 

7 What is needed to form cusps 

Do singular cusps that appear in numerical experiments actually form in 
Nature? To answer this question, it is necessary to understand the structure 
of cusps and what their formation conditions should be. 

As we discerned from the isothermal sphere example, DM particles freely 
intersect the core region r < tq and do not stay there. This conclusion is valid 
for much more numerous halos with a small central mass {a < 1), because 
the gravitational force tends to zero at small radii: 

GM(r) 1 „ 

— oc r^-° ^0 as r ^ . (20) 

Conversely, more massive cusps {a > 1) gravitationally influence DM parti- 
cles and keep them inside the halo. Clearly, velocities of the kept particles 
must be sufficiently low for the particles to remain in the cusp region. There- 
fore, to form a cusp, it is necessary to collect many cold particles in its center. 
At least two questions arise. 

• Where can cold particles in the perturbation field be found? 
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Figure 7: (See color version at www.ufn.ru). The mean halo density profiles 
in numerical models [51] • Dwarf galaxies are shown in red, intermediate 
mass galaxies are shown in green, and galaxy clusters are shown in blue. 
Numbers show the power-law slope exponents din p/dlnr; the normalization 
corresponds to the slope point —2. 

• How can cold particles be transported to the central part of the halo 
without being heated? 

Because the fraction of energy removed during halo formation does not 
exceed 10-15% and the degree of matter contraction during relaxation is 
finite, both these questions should be first of all addressed by numerical 
programs for N-body calculations: where is the 'cooler' of DM particles and 
why is their 'warming up' by small-scale cosmological motions of matter not 
taken into account in calculations? 

We know that in the absence of small-scale perturbations, cusps are 
formed in models of quasispherical collapse (see Section 6). In these an- 
alytic calculations, all particles are initially focused toward the center (by 
the chosen initial conditions), and a significant fraction of the particles, by 
crossing the center, cools due to violent relaxation and remains inside the 
cusp. But in real cosmological conditions, DM particles are 'defocused' by 
small-scale perturbations and hence do not pass through the center. In ad- 
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Figure 8: Demonstration of the cusp problem in galactic DM lialos (data 
from [55]). The symbols show the observed galactic density profiles in (a) 
DDO 47 and (b) ESO 287-G13, the curves show the mean density profiles 
obtained in numerical modeling. 

dition, the fall of DM clumps and their merging 'stir up' the central region 
and hamper effective particle cooling. 

Of course, the quantitative characteristics of these mutually opposite pro- 
cesses must be calculated and checked numerically; however, modern numer- 
ical modeling is not perfect and needs to be improved. 

First of all, it is required to control DM motions in a wide range of scales, 
from large cosmological scales on which the initial particle distributions are 
given, to very small ones, which allows the inner regions of the forming halos 
to be resolved. In modern numerical models, the mass resolution, as a rule, 
does not exceed 10'' Mq, which is insufficient to describe the structure of DM 
cores. Due to the restricted resolution in numerical models, the initial parti- 
cle distributions have no small-scale perturbations, i.e., there is an effective 
cut-off of the perturbation spectrum at short wavelengths. The consecutive 
resolution of selected volumes by the addition of 'new points' inside them 
at the late stage of the evolution (the so-called adaptive method) does not 
help either, because this procedure does not add necessary information on 
the initial small-scale field of cosmological perturbations. 

A more sophisticated method is sometimes used that involves consecutive 
recalculation of the model upon adding new particles in the future dense halo 
regions. This procedure enables small-scale perturbations to be included in 
the calculation in restricted volumes. But such a recalculation complicates 
the problem and strongly increases the computation time. The problem of 
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cusps has not yet been solved. Perhaps this example shows a limitation of 
the CSM, in which DM particles are considered to be cold. In the CSM 
framework, we cannot resolve the contradiction between observational data 
and the results of numerical calculations. The point is that cusps emerge 
during relaxation of the initially cold matter, in which large-scale particle 
flows self-intersect to form caustics, from which a cusp is later formed. But 
the real DM can be initially warm with small random particle velocities. In 
that case, caustics smear out, and cores with a smoother density distribution 
are formed in halo centers instead of cusps. 

Warm DM can be regarded as an additional parameter of the cosmological 
model. The mass of warm DM particles that is required to solve the problem 
of cusps must be not large (~ 10 keV), which allows the particles to preserve 
the residual thermal velocities associated with the hot phase of evolution 
in the early universe. Another possibility of the initial 'heating up' of DM 
particles can be related to the presence of excess power in the perturbation 
spectrum at small wavelengths. We recall in this connection that the spec- 
trum is fairly well determined on scales above several Mpc from observations 
of the cosmic microwave background and large-scale galaxy distribution. At 
the same time, at small scales, the spectrum of initial perturbations can 
deviate from a power law, which can in turn affect the structure of small 
galaxies. 

These assumptions, however, go beyond the CSM framework. Our ap- 
proach is different. We believe that there are currently no strong grounds 
to complicate the fundamental model, and the problem can have a simpler 
solution, which does not require modification of the CSM. The influence of 
small-scale velocity and density perturbations on the internal halo structure 
can be studied analytically using entropy function (8). The use of entropy 
instead of density is motivated by the fact that for most particles, the en- 
tropy increases during the halo formation and relaxation and is an integral 
characteristic of the entire halo formation history. Evaporation of particles 
during relaxation decreases the entropy, but this effect is rather small. 

By definition, DM particles are initially cold, i.e., their entropy is zero. 
On the other hand, there are random small-scale flows and matter clumps, 
which are tidally destructed and scattered when approaching each other in- 
side collapsing protohalos to form folds and intersections, thus heating DM 
particles and chaoticizing their velocities before the violent relaxation begins. 
Therefore, during the first contraction, the small-scale velocity perturbations 
in the protohalo are transformed into chaotic motions of DM particles. This 
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dynamical time 
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time of collapse 



Figure 9: Motion of particles in (a) an ideal gas and (b) a DM protohalo [2]. 
The entropy of particles in the region of size r is determined on the dynamical 
time scale coincident with the collapse time: r/f ^ H^^. 

process can be described in terms of coarse-grained entropy, in analogy with 
the description of the ideal gas [41] . 

8 Ensemble of protohalos 

To introduce the entropy characteristic of particles in a protohalo, the notion 
of an ensemble is needed. It can be introduced as a subsample of the field of 
linear perturbations embracing all spatially bound regions that collapse by 
the time corresponding to the current redshift z. A schematic comparison 
with an ideal gas is shown in Fig. 9. 

A protohalo is characterized by the mass M/j first collapsing at the red- 
shift z, which is called the instant of halo formation. The mass Mji corre- 
sponds to the linear comoving size R: 



A sphere with the radius R = 10 Mpc comprises the typical galaxy cluster 
mass Mio ^ 1.6 ■ 10^^ Mq (see Section 3). 

To single out a protohalo from the common density field, it is necessary 
that the protohalo have an initial threshold mass excess in the sphere of 




(21) 
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radius R which provides the total negative binding energy of this region and 
ensures its collapse by the given time instant (see Appendix A): 

Sr 5c (z) 

iyR = — >u = . (22) 

The function 5c{z), which can be calculated analytically, is controlled numer- 
ically. In the CSM, 5c — g{z) • 5c{z) — const ~ 1.67, where g{z) is the growth 
factor of linear density perturbations (^(0) = 1). The initial ensemble of 
protohalos of interest here, with initial masses 

Mh > Mr (23) 

includes only those space points x for which the mean matter density inside 
a sphere of the radius R around the given point increases the threshold value 
in (22), where 5{x) = p(x)/pni — 1 is the linear density contrast relative to 
the mean density in the Universe p^', the smoothed density contrast 5r with 
the filter of radius R is 



6r = 6r{^) = J (5(x') Wr{\^ - x'l) rfx' , 



(24) 



where WR{r) is the 'window' function 

. . 3 f 1 , r<R 
WR{r) 



47ri?3 [ , r> R 

Clearly, function (24) contains no information on perturbations on scales 
smaller than R. The variance of the smoothed density linear contrast aR, 
linearly approximated to the present moment z — 0, is related to the power 
spectrum of inhomogeneities P{k) by the integral relation (see Appendix A) 



JR 



/•oo 

(5^)^ / P{k)W''{kR)edk, (25) 
Jo 



where W{y) = 3y^'''{smy — ycosy) is the Fourier transform of the window 
function (y = kR). It is known from observations that au ~ 0.8 for the 
radius i? = 11 Mpc, comprising the galaxy cluster mass Mn ~ 2 • lO^'^ Mq. 
In astronomy, it is equivalent to the historically established normalization 
to 'sigma eight', where the index 'eight' corresponds to the physical radius 
R = 8/0.7 Mpc ~ 11 Mpc. 
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We use the obtained ensembles to describe the mean 'thermodynamic' 
characteristics of DM in protohalos. The total entropy of a virializcd halo 
integrates the joint effect of all nonequilibrium processes that occurred dur- 
ing the entire halo formation time, including both the background particle 
(coarse-grained) entropy related to the initial cosmological perturbations and 
the acquired entropy generated during the violent relaxation and hierarchical 
clustering of matter. The cosmological part can be calculated analytically, 
and to estimate the nonlinear component, we use numerical results. 



9 Prom density profile to entropy 

First of all, it should be understood how the entropy is distributed inside an 
equilibrium halo. To simplify calculations, we restrict ourselves by consider- 
ing the power-law approximation of the central halo density profile (5): 

p{r) oc r"" , (26) 

M{r) ^An p{x) x^dx oc r^~" , 
Jo 

with the slope exponent a G (0; 2.5). In numerical calculations of the halo 
formation, profiles with a > I emerge, while observations typically show that 
a < 1 (Fig. 10). In the first case, we are dealing with a cusp (1 < ck < 2.5) 
and in the second case, with a core (0 < a < 1). 

The physical difference between the cusp and the core follows from the 
behavior of the effective pressure, which can be determined from Eqn (7): 

P(r) = ci + C2r2(^-"), (27) 

where Ci and C2 are integration constants. The critical slope is a = 1: in a 
core, the pressure is finite, while in a cusp it diverges at the center. 
From Eqn (8), we obtain the entropy mass function 

E(M) oc ci + C2 oc , (28) 

P 5a ^ 6 — a 

3(3 — a) ' 3(3 — a) ' 

where the parameter /? = rflnE/rflnM lies in the range between /3i and ^2 
and its local value depends on the current mass of the halo. At a — ac= 1, 
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a 

Figure 10: Distribution of the power-law exponent a in the central part of 
galaxies with a low surface brightness (see [61]). 

the interval of possible values of /3 contracts to the point 

/3l = /32 = /3e = ^. 



We have passed from the radius r to the comoving mass M in the sphere of 
radius r because this mass is conserved during both linear perturbations and 
the halo relaxation. 

As stressed above, the entropy mass function accumulates the action of 
irreversible processes during DM evolution and determines the equilibrium 
halo density profile. In the considered range of a, we have /3i_2 > 0, and 
hence E{M) — )■ as M — )■ (Fig. 11). But in the cusp center, the entropy 
of particles is very small, 

/3>5/6: M-^/^E(M)-^0 as M^O, 

while in the core center it is much higher, which makes the cusp formation 
impossible, 

/3 < 5/6 : E(M) ^ oo as M -> . 

For the critical value /3 = 5/6, it follows from (8) and (28) that 

M oca\ (29) 

This formula reproduces the well-known empirical Faber-Jackson |62] and 
TuUy-Fisher [53] relations, which relate the observed velocity dispersion of 
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Figure 11: Profiles of (a) tfie density and (b) tlie entropy in a DM lialo [2]. 

stars or gas to tlie mass (luminosity) in elliptical galaxies or spherical sub- 
systems of spiral galaxies. 

Figure 12 illustrates the difference between the cusp and the core of an 
equilibrium halo from the standpoint of the entropy distribution. It follows 
that the formation of the central cusp or core depends on the total entropy 
of matter acquired before the formation of the halo, during and after its 
formation. 

We can now quantitatively answer the question as to how the initially 
cold DM particles of the CSM acquire a nonzero temperature. 

The halo formation process can be divided into several consecutive stages. 

• Initial DM motions are potential and are correlated, relative veloci- 
ties of close particles are small, and the DM velocity field is smooth. 
Random deviations from the mean velocity determined by the initial 
perturbations spectrum increase with increasing the scale to saturate 
at relatively large scales (~ 40 Mpc) (see Fig. 14 in Section 10). 

• During contraction of a protohalo, particles approach each other and 
mix at small scales, and the initial phase of the small-scale part of the 
random velocity is lost: DM particles heat up (the background entropy 
mass function Eb arises) H 

• The chaotization of the regular contraction rate and transformation of 
the kinetic energy of the contraction into thermal energy occur some- 

Coarse-grained entropy can be defined as the fraction of the phase space volume per 
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M 



Figure 12: The entropy mass function of an equilibrium halo. Curve 1 shows 
the cusp, curve 2 shows the core, and the dashed curve corresponds to (3c = 



what later during violent relaxation and merger of clumps (the entropy 
Eg is generated). 

• The resulting density profile of equilibrium halos is obtained by adding 
the initial entropy Ef, and the acquired entropy Eg. 

In Sections 10 and 11, we obtain the initial mean entropy profile (E{,) using 
the corresponding dispersion of linear velocity perturbations [see Section 11, 
formula (49)]. This representation assumes an almost adiabatic contraction 
of matter and gives the lower limit of the mean function Eb{M) at small M. 



where AiV is the number of particles in the volume Ax = Axi Aa;2Aa;3, including the 
oppositely directed particle flows, Avi is the relative particle flow velocity on the scale 
Axi (Av = A?;iAu2Au3), and uav is the dispersion of relative particle flow velocities 
inside the volume Ax. The coarse-grained entropy does not decrease during relaxation 
of particle flows and is transformed into the real entropy mass function of the virialized 
system. 



5/6. 



particle: 
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10 Initial cosmological perturbations 

The field of scalar cosmological perturbations is characterized by three gauge- 
invariant variables: the space vector of the shift of a material point from the 
unperturbed position 

S = S(x) , (30) 

the total velocity of matter 

V = V(^,x), (31) 
and the comoving density perturbation 

^ - 1 = giz) . 5(x) . (32) 

The dependence of the mean density on the redshift has the form p{z) = 
Pm (1 + 2;)^, where the DM density is given by (see Section 3) 

ottG sm-^ 

The mean particle shift in the Universe is zero, and its dispersion is deter- 
mined by the integral of the power spectrum of density perturbations: 

(poo \ 1/2 

J P{k)dkj ~14Mpc. (33) 

A sphere with the radius R = as contains the mass of a rich galaxy cluster 

Ms = — p^alc^ A- 10'^ Mq. 

Relations between variables are determined by the linear perturbation 
theory (see [Ij for more details): 

r(^,x) = (l + ^)-^(x + ^S), (34) 

\ = x = H{x- g'^S) , 5(x) = - div S , (35) 



y-Hx=-Hg:,S, av = /P)^ 200 J^^, (36) 
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where t{z, x) and x are the Eulerian and Lagrangian coordinates of the 

medium, v = Vpcc and dv are the pecuhar velocity and its variance, and 
is the derivative with respect to the redshift z. The growth factor g = g{z) 
is normahzed to unity at zero redshift {g{0) = 1), and at 2; > 1, 

1.3 

9 - 



1 + z 

The Hubble function and the relation between the redshift and the peculiar 
velocity of matter at ^ > 1 have the form 

H c^0.5-{l + zfHo, (37) 



r. rr I o 45 / S \ km 

V ~ 0.5 • E^^^Y^zg^ ~ ^== t:^ — 



To characterize small-scale motions of matter in a protohalo, we introduce 
the vector of the relative shift of medium points separated by a distance r 
[see (30)]: 

, , S(x + r)-S(x) 

s - s(r, x) = ^ -1 ^ . (38) 

V2crs 

The vector s describes the field of the normalized peculiar velocity of matter 
with the current radius vector r referenced to the point x. We are interested 
in the mean relative velocity and its variance during averaging over different 
ensembles (or points x) in different regions of the Universe. 

The modulus of the radius vector r — |r| is fixed in the interval of linear 
scales corresponding to inner regions of a gravitationally bound halo of size 
R: 

r<R. (39) 

At R < as, these scales relate to the short- wavelength part of the spectrum 
P{k) (to the right of its maximum in Fig. 13), 

/co-s > 1 : P{k) oc k-^hi^{kas) ■ 

Hence, we find the estimate for ar [see (25)] 

POO 

(jI= / P{k)W'^{kr)k'^dk(xhi^{2(7s/r). (40) 

JQ 
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Figure 13: The power spectrum P{k) in the CSM, crio = 1. 

Taking the power spectrum normahzation a^o ~ 1 into account, we can 
approximately calculate 

1 < (7fl ~ ln3/2(2(7s/r) at R < 10 Mpc . 

At > 10 Mpc, a more accurate estimate for au < 1 should be used, taking 
the change in the spectral slope P{k) at scales of rich galaxy clusters into 
account. 

Averaging over the entire Universe, we obtain the mean relative velocity 
and its variance: 

(s) = 0, (41) 

- is') - I^^Pik) (l - ^) dk . (42) 

The random part of the velocity as = Os{r) monotonically increases with 
increasing the scale (Fig. 14). At r < us, its growth is of a quasi-power form 
with a smoothly decreasing slope exponent, starting from unity: 

ra,. 



V6 



The growth in the function cXs decreases at r > cxs and saturates at r ^ cxs: 
as^l. At r = as, ~ 1/2. 
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Figure 14: The variance of the relative velocity o"g(r) and its longitudinal aj^ 
and transverse a\ components (see j2]). 

The correlation radius of function (42), at which crg(rs) = 1/2, is ~ 
36 Mpc. We note that this is the mean correlation radius of the relative 
velocity of matter. The variance of vector (38) is anisotropic (different for 
the longitudinal and transverse components) and depends on the location of 
particles relative to the halo center r = re : 



where the tensor pij = Sij — CiCj is the projector on the plane perpendicular 
to the vector e. 



The correlation radii of the longitudinal and transverse variance are 20 Mpc 
{al - 1/6) and 49 Mpc {a^ - 1/6). 

Averaging over the ensemble of protohalos (see Appendix B) with a den- 
sity contrast [see (22)], we obtain the conditional expectations of small- 
scale velocities of matter in a protohalo: 



= = Or|| + 2(T_L , 



(43) 



(44) 





(s| Vr) = CZ/r, 



(45) 
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a%^{s'\u^)-{s\u^Y = al-c\ (46) 
where the cross-correlation coefficient has the form [see (39)] 

c={sun) = 7=^ I P{k)W{kR)W{kr) Pdk - "^^^ 



3^2 o"s ctr Jo 3v^ 



Function (45) determines the universal velocity proffie of a collapsing proto- 
halo with the mass M^. Conditional variance (46) describes dispersions of 
random deviations of the velocity from the mean value inside the protohalo: 



11 A solution to the problem of galactic cusps 

The suppression factor of peculiar motions of DM in a protohalo relative to 
the motions of matter in the Universe whole, 

f[r\R) = ^ ~ 1 - — 

is shown in Fig. 15 for the galactic mass 1.3 ■ 10^^ Mq {R = 2 Mpc). For 
galaxies and groups of galaxies, the cross-correlation inside the protohalo is 
insignificant: 

^s%(0 ^ ^ (^^r - ^ 4 , (48) 

and the variance of small-scale velocities of DM has the universal form, which 
is independent of the protohalo mass [see (36), (37)], 

however, the value of the variance depends on the halo formation instant z. 

By assuming that the coarse-grained entropy E{M) does not change in the 
protohalo contraction (see footnote 3 in Section 9), we obtain the mean profile 
of the background entropy mass function (8) in parametric representation [21 

2 / \ 2/3 



2 



{E,) = ^{-^] oca,^(r), M = ^ r\ (49) 
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Figure 15: The suppression factor f{r\R) of peculiar velocities inside a halo 
with the mass M = 1.3 x 10^^ Mq (linear size R = 2 Mpc) (see |2]). 



Table 1: The power-law exponent of the background entropy function /3b for 
the central halo masses M = 10"M„ < Mr? < 4 ■ lO^^M^. 



n 


12 


10 


8 


6 


4 




0.34 


0.45 


0.50 


0.53 


0.56 



Hence, we calculate the power-law slope exponent /3f, depending on the index 
n of the central halo mass M = W Mq < Mr < Ms ^ 4: ■ 10^"^ Mq : 

,_ dHT^,) _2 1 ^2 1.3 

dlnM 3 ln(2(Ts/r) ~ 3 16 - n ' ^ ' 

The Table lists numerical values of Pb for different n. We see that in the 
halos of galaxies and groups of galaxies, the background entropy hampers 
the central cusp formation, because 

/Sfc < 0.6 < ^ . (51) 
6 

As R increases in the region R> as and the inner radius r < R increases pro- 
portionally (i.e., for large halos with Mji > 4- lO^^M©) the cross-correlation 



We note that the radius of influence of dark matter r cannot be considered small, 
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(s vr) cannot be neglected anymore [see (47)], and small-scale motions of DM 
in the protohalo are suppressed. Therefore, the effect of the background en- 
tropy of DM inside rich galaxy clusters decreases. 

Thus, DM cores are formed in the central region of halos of galaxies and 
galaxy groups. Following astronomical tradition, these DM cores can be 
appropriately called a spherical subsystem or 'bulge' of DM (in analogy with 
the baryonic bulge). By assuming that a ~ 1 in the region with maximum 
luminosity [see (29)], the Faber-Jackson relation must be vahd for bulges, 
which is actually observed in elliptical galaxies. In disk and spiral galaxies, 
the bulge occupies a relatively small part of the halo, although its luminosity 
in the Ha line can be significant, which apparently explains the TuUy-Fisher 
relation. 

In large spiral galaxies, the DM bulge is deformed due to baryonic inflow[f] 
however, its size is still quite large. For example, the bulge of our Galaxy 
occupies several kiloparsecs. As r increases, it transforms into a DM distri- 
bution with p oc and an almost flat rotation curve, which in turn extends 
up to 10 kpc and, possibly, even further beyond the disk edge; further DM 
distribution can be probed by the motion of dwarf satellites. 

On the other hand, numerical modeling shows that in gravitationally 
bound galaxy groups, the central bulge could be formed only in compact 
groups. There is a large class of smoothed galaxy groups that are only 
partially relaxed; such is the Local Group of galaxies. Its mass ~ 2 ■ 10^^ Mq 
is mainly due to the masses of two large spiral galaxies (the Milky Way and 
Andromeda Nebula) and also includes a small addition from masses of several 

taking the cooling baryon flows and massive cD-galaxies in the cluster centers into account. 
For a similar reason, it makes no sense to consider central DM masses lower than that 
of black holes in galaxy centers. We should bear in mind in this connection that the 
characteristic scales of DM distributions are much larger than the sizes of central baryon 
condensations. 

^It should be noted that during a slow baryonic inflow, the DM cusp (if it exists) 
would not be destroyed, but only enhanced. In this sense, it is asymptotically stable. 
The question is: what is the DM density in the halo center? This density, in particular, 
determines the possible flux of DM annihilation products (photons, positrons, antiprotons, 
etc.). Additional studies are needed to answer this question. Dwarf galaxies with low 
surface brightness appear to be ideal objects for studies of the central DM distribution 
(see Fig. 3), in which star formation is suppressed and the gas sustained in rotational 
equilibrium does not fall onto the halo center. Studies of spiral galaxies have not provided 
evidence of the existence of central cusps, either. Weak gravitational lensing plays an 
increasingly important role in reconstructing the DM surface density. 
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dozen dwarf galaxies. 

When comparing theoretical predictions with observations, a wide dis- 
persion in the background entropy distribution (and the corresponding inner 
density profiles) of relaxed halos, which is actually observed in galaxies (see 
Fig. 10) and follows from the CSM theory, should be taken into account. 
Indeed, for the initial Gaussian velocity perturbations, the probability dis- 
tribution of the background entropy has the form 

«A) = !5fPt4Z2)rf/„ (52) 

where 

We see that the halo background entropy has a wide distribution and its 
variances around mean value (49) are large: 

ih) = 1 , im = 3(A) = 3 . 

These large values 

varEf, > (Eb)^ 

imply significant deviations from the mean power-law slope exponents of 
the central density profiles in different galaxy halos, which is confirmed by 
observations. 

The problem of galactic rotation curves is tightly connected with the inner 
structure of equilibrium halos, and we discuss it in Section 12. 



12 Galactic rotation curves 

In Section 5, we considered fiat rotation curves of stellar and gas disks ob- 
served in many galaxies, which suggest an isothermal DM distribution inside 
halos in the rage of radii and masses <10. The corresponding power- law ex- 
ponents of the density and entropy profiles in such systems [r > Tq; see (11) 

and (24)] are 



a ~ 2 , 




As noted in Section 6, similar distributions are generated in analytic models 
of the spherically symmetric or quasispherical collapse of gravitating dust 
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matter with a smooth initial density distribution, in which background small- 
scale perturbations are definitely absent. During nonlinear contraction of such 
a cloud, caustics (self-intersections) and multi-fiow motions of matter arise in 
its central region, inducing a violent, coUisionless relaxation. As a result of 
the gravitational redistribution of energy between spherical shells of matter, 
power-law density and entropy profiles are formed in the central parts of the 
halo, with the exponents 

1.6<ag<2, l<l3g<1.3. (53) 

It follows from (53) that the violent relaxation leads to the entropy dis- 
tribution implying the formation of a cusp. However, these exponents are 
related to rather large radii. In the central halo, generated entropy (53) can- 
not compete with the background one and turns out to be negligibly small 
(see the Table). Therefore, the background entropy effect constrains the 
central density and hinders the formation of cusp-like halo profiles. 

The same conclusion is valid for the numerical experiment modeling of 
the most complex stages of halo formation that cannot be treated analyti- 
cally. Violent relaxation, in which radial trajectories play the dominant role, 
is not the only and, possibly, the main relaxation process of coUisionless sys- 
tems that decoupled from the Friedman fiow during cosmological expansion. 
An important role here is played by anisotropic collapse and the effects of 
merger of clumps with different masses that were formed before the central 
massive part and simultaneously with it. All these processes, which are actu- 
ally observed in nonlinearly relaxing gravitationally bound systems at high 
redshifts, are called hierarchical clustering. They include the tidal binary 
(and also triple, etc.) interaction of clumps of matter captured in the com- 
mon gravitational field and their scattering and merging into more massive 
objects with the subsequent repetition of the entire cycle. 

The development of numerical modeling techniques for a large number 
of gravitationally interacting bodies has led to the notion of the universal 
halo density NFW profile, which appears after averaging the density profiles 
in several hundred clumps of different masses and sizes formed during the 
merger. The initial particle distribution in the N-body problem is taken 
to be close to cosmological conditions, but because the number N is finite, 
the dynamical range of model scales is restricted. As noted in Section 6, 
the NFW profile asymptotically approaches a cusp-like power law with the 
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exponents 

5 

Q^g = 1 , /3g = g • 

In numerical modeling of the iV-body problem, the background entropy of 
small-scale perturbations is apparently underestimated due to the restricted 
spatial resolution of this method. Hopefully, as the power of numerical meth- 
ods increases, the background entropy effects can be properly taken into 
account. 

Already much evidence has appeared showing that with increasing A^, 
the exponent a decreases below unity; however, these are single examples 
and no statistics can be applied. Due to the small-scale cut-off of the model 
power spectrum, the initial entropy effect is suppressed in comparison with 
the analytic estimates shown in the Table. 

Of course, the background entropy effect on the internal halo structure 
and the formation of the central core with finite pressure and DM density 
should be confirmed by numerical experiment. Nevertheless, today we can 
already verify the GSM predictions by measuring velocities of gas and stars 
in galaxies and galaxy clusters. 

Figure 3 demonstrates galactic rotation curves that cannot be described 
by formula (19) but are well approximated by a smoother profile in (18). 
Such examples are typical for galaxies (see, e.g., |6^). In the central region 
with r < r^, the fit constructed using the observational points is similar to 
the isothermal sphere approximation [see (16) and (17)]; but at large radii, 
the density decrease is inversely proportional to and circular velocities 
slowly decrease. The circular velocity reaches a maximum at the point of the 
local density profile with the exponent a = 2. 

Simulated NFW profiles meet with another problem: they do not repro- 
duce fiat galactic rotation curves in the observed range of scales. Neverthe- 
less, numerical profiles provide a good approximation of the mass distribution 
in galaxy clusters (see Fig. 4, 5). 
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13 Effect of the background entropy 



To understand the physical meaning and variety of the observed rotation 
curves, the initial background entropy of DM in nonlinear halos should be 
taken into account. 

To clarify the physics of this effect, we consider the simplest model of the 
joint entropy mass function in the form 



where the parameter f{, e (0, 1) expresses the relative contribution of the 
background entropy at a fixed inner halo radius corresponding to the circular 
velocity maximum or the local density profile slope a = 2. The functions 
Eft_g(M) have a power- law form with the corresponding slope exponents in 
the ranges 



We now can calculate the rotation curves of particles in equilibrium DM 
halos analytically and compare them with observations and numerical model 
calculations. The results are presented in Fig. 16. The rotation velocities 
Vrot = v(r) as functions of the inner radius r are obtained for models of hier- 
archical {^g — 5/6, Fig. 16a) and violent {^g — 4/3, Fig. 16b) relaxations for 
two limit values of f;,, f;, <C 1 and f^ ~ 1, and three values of the background 
entropy indices. Velocities and radii are normalized to the corresponding 
maximum values of the velocity and the radius r_2 at which it is attained: 
Vmax = v(r_2) > v(r). We see from this figure that for halos with (5^, < 0.5 
and any parameters and f^ in the chosen intervals, the induced rotation 
curves fully cover the range between the NFW and Burkert approximations. 
For fb <^ 1, the curves shift to the NFW region, while as f;, increases, they 
become similar to the Burkert profiles. 

It follows that the observed rotation curves are successfully reproduced 
by simple models (54) with the corresponding fitting parameters f^ and 
Distributions of the background (52) and generated (55) entropies in the 
halo allow explaining the required diversity of the observed galactic rotation 
curves. As the exponent (5^ increases, the rotation curves concentrate closer 
to the NFW profiles. 

The examples considered here prove the following. 





(54) 




(55) 



37 




Figure 16: Normalized rotation curves v{r) for models of (a) hierarchical 
{f3g = 5/6) and (b) violent {f3g = 4/3) relaxation at f;, ^ 1 (the dashed curve) 
and ffe ~ 1 (the solid curve) (see |3I])- The NFW and Burkert approximations 
are respectively shown by X's and the dotted line. 

• The background entropy prevents the formation of central DM cusps 
in halos in the mass range 10^ — 10^^ Mq. For heavier and lighter halos, 
the background entropy effect is attenuated. 

• Taking the effect of the background entropy on the density distribution 
inside a halo into account allows reproducing the rotation curves in a 
wide range of scales and can solve the problem of DM cusps in the CSM 
framework. 

To conclude, we note that the models considered here are based on simple 
assumptions of the DM particle distribution in the halo and neglect the 
contribution of baryonic matter. The full problem of the evolution of galactic 
density profiles, including DM and baryons and accounting for the effect 
of different dissipative processes on the angular momentum transfer from 
the central parts of galaxies, is undoubtedly very complicated and requires 
further studies, both analytic and numerical. 
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14 Conclusion 



The entropy method of description of viriahzed DM systems considered here 
allows an analytic treatment of these complex nonlinear structures. It enables 
connecting the inner DM density profiles with characteristics of both the ini- 
tial small-scale field of density perturbations and the nonlinear large-scale 
relaxation of gravitationally contracted matter. We conclude that cosmo- 
logical random motions of matter 'warm up' DM particles in the collapsing 
protohalos. We have shown that taking this effect into account 

• leads to suppression of cusp-like density profiles in the forming halo 
and to the formation of DM cores in galaxies; 

• allows explaining the diversity of galactic rotation curves, both ob- 
served and obtained in numerical simulations; 

• helps to solve the problem of the inner structure of equilibrium DM 
halos in the CSM framework. 

The obtained analytic results should be confirmed by numerical iV-body 
experiments, which can be done in the future with improved spatial resolution 
of the central halo parts. 

The authors thank S V Pilipenko for the fruitful discussions. 
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15 Appendices 

Appendix A. Correlation functions of the initial pertur- 
bations 

Fourier harmonics of the linear density perturbation field are 5-correlated: 

5(x) = (27r) -3/2 J 4 e^^'dk , (5k5k') = ^Tr^Pik) 6^^^ (k - k') , 

where the averaging (...) is done over all realizations of the random Gaussian 
field, 5^3) (k) is the three-dimensional Dirac 5-function, and P — P{k) is the 
power spectrum. 
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Using this decomposition, we can construct any correlators of the cosmo- 
logical perturbation field, for example 



(5(x + r)5(x)) 



(S(x + r)5(x)) 



-- r 



^ , , sin kr ^ 

Pik)— k^dk, 

kr 

P{k)W{kr) k'^dk, 



(,S,(x + r)5,(x)) = e,efe / P{k) 

Jo 



sin kr 2 

1 



dk + 



+ oPjk / P{k)W{kr)dk 
■J Jo 



where e = ej = rj/r is the unit vector in the direction r, pjk is the projection 
tensor [see (43)], 6 — — divS, and S is the shift vector of matter elements 
[see (30)-(37)], 

We can represent smoothed fields similarly, for example, 
Sr = {2tt)-^/^J 4 W{kR) e^^dk , 

(S(x + r) 5k(x)) = / P{k) W{kR) W{kr) k^dk , 

3 Jo 

etc. 



Appendix B. Conditional probability distribution 

To obtain the conditional probability p(s|z/i?) of s for a given value of z/r, 
we consider the Gaussian distribution of two variables xa = (s, i^r) with the 
index A ranging over the values i, R: 

. . exp( — 2 AR 

P ^ ,2 /77 ' K = C XaXb , 

{Zn) V C 



c 1 



where c^^ is the inverse matrix to the matrix [see (43), (45)] 

Cab = {xaXb) = 
40 



C = det {cab) = c^l , (7|| = (T^ - . 
After simple transformations, we obtain 

— k'^ + = c^^SiSj , s = s — ciyR , 

where c'-^ is the matrix inverse to 

Cij = {siSj\pR) = a-|| ejCj + a]_Pij . 

Hence, we find the conditional probabihty distribution of the vector s and 
the variance of velocities inside a halo with the density contrast i/r: 

^ _ Pj^^ ''r) _ exp(-ft:V2) 

2 /-2| \ 2 2 
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